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ABSTRACT 

The study of elemental abundances in Damped Lyman Alpha systems (DLAs) at high 
redshift represents one of our best opportunities to probe galaxy formation and chem- 
ical evolution at early times. By coupling measurements made in high z DLAs with 
our knowledge of abundances determined locally and with nucleosynthetic models, we 
can start to piece together the star formation histories of these galaxies. Here, we dis- 
cuss the clues to galactic chemical evolution that may be gleaned from studying the 
abundance of Co in DLAs. We present high resolution echelle spectra of two QSOs, 
Q2206-T99 and Q1223+17, both already known to exhibit intervening damped sys- 
tems. These observations have resulted in the first ever detection of Co at high redshift, 
associated with the z a t, s — 1.92 DLA in the sightline towards Q2206— 199. We find that 
the abundance of Co is approximately \Z Q and that there is a clear overabundance 
relative to iron, [Co/Fe] = +0.31 ± 0.05. From the abundance of Zn, we determine that 
this is a relatively metal-rich DLA, with a metallicity approximately \Z Q . Therefore, 
this first detection of Co is similar to the marked overabundance relative to Fe seen 
in Galactic bulge and thick disk stars. 
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1 INTRODUCTION 

Our understanding of the processes associated with galactic 
chemical enrichment, such as stellar nucleosynthesis, dust 
depletion and mixing mechanisms, can benefit greatly from 
the comparison of local abundances with those measured at 
high redshift. At low z, abundances can be determined for 
the Milky Way and its near neighbours from stars, H II re- 
gions and the interstellar medium (e.g. Chen et al. 1999; 
Kobulnicky and Zaritsky 1999; Shetrone et al. 2000; Izo- 
tov and Thuan 1999; Savage and Sembach 1996). At high 
redshift, the most accurate and extensive abundance mea- 
surements are those determined from Damped Lyman Alpha 
systems (DLAs). These systems are the highest H I column 
density members of the quasar absorption line menagerie 
and are traditionally defined as having N(H I) > 2 x 10 20 
atoms cm~ 2 . Coupling measurements of metal abundances 
in high redshift DLAs with our knowledge of local enrich- 
ment patterns plays a crucial role in our understanding of 
the star formation histories of these distant objects and the 
early stages in galaxy evolution. 

In the context of investigating galactic chemical evo- 
lution, several key elements have previously been targeted 
as particularly informative diagnostics. The utility of these 
key metals is often based on our knowledge of their relative 



production and release timescales in stars (and subsequent 
supernovae) of different masses. For example, the a elements 
(Mg, Si, Ca etc.) have been extensively used to investigate 
the amount and rate of massive star formation. An over- 
abundance of a relative to Fe in metal-poor Galactic halo 
stars chronicles the epoch of early star formation (Tinsley 
1979; Ryan, Norris and Beers 1996) and in the bulge the 
same overabundances in solar metallicity stars are indica- 
tive of a 'get rich quick' population (McWilliam and Rich 
1994; Rich and McWilliam 2000) i.e. one which is as old as 
the Galactic halo but as metal-rich as the Galactic disk. In 
DLAs, there is on-going discussion about the abundances of 
a elements, for which nucleosynthetic origins of abundance 
patterns must be disentangled from dust depletion effects. 
Although there is some evidence that Si may be overabun- 
dant in DLAs (Lu, Sargent and Barlow 1998), [Si/Zn] is solar 
in at least some systems (Pettini et al 2000) . Moreover, since 
Si is highly refractory, the interpretation of [Si/Zn] measure- 
ments relies on a dust correction factor based on [Zn/Cr], 
which is quite uncertain. Conversely, the abundance of S, an 
a element not depleted onto dust, has been found by Cen- 
turion et al. (2000) to be underabundant relative to Fe, a 
result at odds with the Type II supernova pattern observed 
in the Milky Way. These results have been interpreted as 
possible evidence for star formation having occurred at a 
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low rate or in sporadic bursts. However, the Centurion et 
al. (2000) data were obtained with resolution R = 5000 and 
need to be followed up with higher resolution data in order 
reliabley distinguish the S lines from the Lyman a forest. 
Further clues to the star formation history of DLAs may 
be gleaned from N/O ratios. The large observed spread in 
N/O in DLAs can be explained by the delayed release of 
primary N (relative to O) in metal-poor intermediate mass 
stars (Pettini et al. 2001; Lattanzio et al. 2001). 

Here we focus on cobalt, an element which although 
relatively widely studied in Galactic stars, has not yet been 
featured in DLA abundance studies. The intriguing trends 
of [Co/Fe] as a function of [Fe/H] in the various Galactic 
stellar populations indicate that cobalt may be another el- 
ement that can provide important clues to chemical evolu- 
tion. In sub-solar metallicity thin disk stars, Co is found 
to be slightly underabundant with respect to Fe ([Co/Fe] 
~ —0.1, Gratton and Sneden 1991), although this may be 
related to their solar log gf scale being on average 0.06±0.04 
higher than the laboratory scale adopted by many other ob- 
servers (see Section 5.1 for more details). This trend con- 
tinues in moderately metal-poor halo stars down to metal- 
licities of [Fe/H] ~ —2.5. For more metal-poor halo stars, 
a considerable overabundance of Co with respect to Fe is 
observed, < [Co/Fe] < 0.8, (McWilliam et al. 1995; Ryan, 
Norris and Beers 1996). Interestingly, this coincides with a 
downturn in relative Cr to Fe abundances, such that Co and 
Cr are closely anti-correlated in the metal poor population. 
Similar [Co/Fe] overabundances are observed in the bulge 
([Co/Fe] ~ +0.3, McWilliam & Rich 1994), which, we re- 
call, is as old as the halo, and to a lesser degree in the thick 
disk ([Co/Fe] ~ +0.1, Prochaska et al. 2000). 

The detection of Co II transitions is a challenging possi- 
bility for ISM absorption line spectroscopy, both locally and 
at high z. Although there are numerous transitions that lie 
at convenient rest wavelengths for both local and high red- 
shift studies, only two interstellar measurements currently 
exist. These detections have been achieved in the local ISM 
in the sightlines towards £ Oph (Federman et al. 1993) and 
p Oph A (Mullman 1998b) using the Co II A1466 transition 
which has an equivalent width in these interstellar clouds of 
only 0.5 and 4.7 mA respectively. One of the advantages of 
moving to higher redshifts is that the relatively strong UV 
lines become shifted into the optical where one has access 
to more efficient detectors and larger telescopes. In addi- 
tion, the observed equivalent width is increased by a factor 
of (z + 1). However, the fundamental problem remains that 
the Co II lines are all intrinsically weak and that the solar 
abundance of Co is 8 x 10~ 8 that of hydrogen. In addition, 
Mullman et al. (1998a,b) have provided new laboratory mea- 
surements for Co II /-values (oscillator strengths), typically 
revising the old values downwards by a factor of approx- 
imately two. This means that the spectral detection limit 
has effectively been increased by the same factor. Finally, 
the two local detections have brought to light an additional 
problem, namely that cobalt appears to be a highly refrac- 
tory element that is depleted onto dust to a similarly high 
degree as iron (Mullman et al. 1998b). 

Despite the technical challenge of observing Co in the 
high redshift ISM of DLAs, it is now feasible with high res- 
olution spectrographs on 8 - 10 m class telescopes to de- 
tect this element, in reasonable exposure times, down to log 



iV(Co) ~ 12.5. Here we present the first ever detection of Co 
in a DLA, at z a b s = 1.92, determined from a spectrum ob- 
tained with the Ultra-violet and Visual Echelle Spectrograph 
(UVES) on UT2 of the VLT. In addition, we present a com- 
bined UVES and HIRES (High Resolution Echelle Spectro- 
graph) spectrum obtained at the VLT and Keck telescopes 
respectively, which can also offer an interesting upper limit 
to [Co/H] in a second DLA. We discuss the implications for 
the star formation history of these high redshift galaxies and 
compare our observations with recent nucleosynthetic yield 
models. 



2 TARGET SELECTION 

In order to at least place interesting limits on the abundance 
of Co in DLAs, it is desirable to reach a sensitivity level of 
[Co/Fe] = 0. Such a detection limit is sufficient at least to 
distinguish between the marked overabundance of Co rela- 
tive to Fe seen in the halo, bulge, and to a lesser degree in 
the thick disk, and the mild underabundance in the Galactic 
thin disk. Therefore, in order to maximise the likelihood of 
a Co II detection, we have selected bright QSOs with fore- 
ground DLAs known to contain high column densities of Fe. 
This criterion can either be achieved in moderate metallicity 
DLAs with high iV(H I) or in lower column density systems 
with relatively high levels of enrichment. The former type of 
system is relatively rare due to the power law distribution of 
column densities in H I absorbers which shows that systems 
with iV(H I) Si 2 xlO 21 cm -2 are very scarce. Similarly, 
few metal rich (Z > jqZq) DLAs have been identified at 
any redshift (Pettini et al. 1999) , indicating that extinction 
due to dust may be biasing current samples (Ellison et al. 
2000). Finally, as dust extinction may be greater in metal- 
rich and high-iV(H) DLAs (e.g. Prantzos and Boissier 2000), 
it is particularly challenging to identify QSOs absorbed by 
such systems that also are bright enough to provide the req- 
uisite S/N. 

Having searched the literature for candidate targets, we 
initially identified two QSOs for study, Q2206-199 (B = 
17.8, z™ = 2.559) and Q1223+17 (B = 18.5, z em = 2.936). 
The former of these targets actually exhibits 3 separate 
DLAs in its absorption spectrum (Prochaska and Wolfe 
1997); here we focus on the system at z a bs = 1.92. Although 
this system only has a moderate column density compared 
with other DLAs (JV(H I) = 4.8 xlO 20 cm" 2 , see section 
4), it is amongst the most metal-rich systems known with 
Z ~ \Zq. In fact, Prochaska and Wolfe (1997) were able 
to detect some elements rarely seen in DLAs in this system, 
such as Ti. Conversely, Q 1223+ 17 is a very high column 
density absorber, N(H I) = 3xl0 21 cm" 2 and therefore has 
a relatively high column density of Fe despite its low metal- 
licity (Pettini et al. 1994; Prochaska et al. 2001). 



3 OBSERVATIONS 

We present high resolution echelle spectra for these two high 
redshift QSOs. The data were obtained with UVES on the 
VLT (Q2206-199 and Q1223+17) and with HIRES on the 
Keck telescope (Q1223+17). A summary of the observations 
is given in Table 1. 
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Table 1. Summary of Observations 



Quasar 


Telescope /Instrument 


Max. A Coverage 


Resolution 


Int Time (s) 


S/N 


Q2206-199 


VLT/UVES 


3300 - 6650 A 


42 000 


17100 


20 - 70 


Q1223+17 


VLT/UVES 


3800 - 10250 A 


45 000 


10 800 


20 - 50 


Q1223+17 


Kcck/HIRES 


4750 - 7200 A 


38 000 


19100 


~ 30 



3.1 Details of Observations and Data Reduction 
of Q2206-199 

Spectra of Q2206— 199 were obtained over 3 nights from May 
28-30 2000 with UVES on the VLT in conditions of good see- 
ing, typically 0.5 - 0.7 arcseconds, although for brief periods 
extremes of 0.4 and 1.3 arcsecs were reached. A slit width of 
1.0 arcsec was fixed throughout the observations, oriented 
along the parallactic angle. The data presented here were 
obtained with the dichroic 1 A cen = 390 + 564 configura- 
tion and read out in 2 x 2 pixel-binned, high-gain mode. 
The blue arm CCD is a single thinned EEV chip whereas 
the red arm consists of a mosaic of a second EEV chip and 
an MIT chip, separated by approximately 0.96 mm. A sin- 
gle dichroic setting therefore has non-contiguous wavelength 
coverage and in the A cen = 390 + 564 setting there are gaps 
between 4525 - 4625 and 5600 - 5680 A. 

The extraction of the spectra was achieved using a 
customized version of the UVES pipeline which is based 
on ECHELLE routines in the data reduction package MI- 
DAS. A detailed description of this process can be found in 
Ballester et al. (2000). The spectra were optimally extracted 
separately for the three CCDs and the variation in the con- 
tinuum flux level, due primarily to the blaze function of the 
echelle grating, was removed. The resolution of the spectra, 
as determined from ThAr arc lines, is R ~ 43 000, or ap- 
proximately 7 km s" 1 FWHM. 1-D spectra were obtained 
by joining the orders and converting the wavelengths to a 
vacuum heliocentric scale. The individual spectra were then 
co-added with a weighting proportional to their S/N. The fi- 
nal step was to normalise the spectrum by dividing through 
by a spline function fitted through absorption-free regions of 
data. Since most of the metal lines which we study here are 
found redward of the Lyman a emission, the continuum can 
be accurately fitted since there are relatively few absorption 
features and the S/N is quite high. 



3.2 Details of Observations and Data Reduction 
of Q1223+17 

The spectrum presented here of the z a i, s = 2.466 DLA to- 
wards Q 1223+ 17 is a composite of HIRES plus UVES data. 
The UVES data were obtained during the same run as the 
Q2206— 199 spectra and reduced in an identical way. The 
only difference for the data of Q1223+17 is that they were 
obtained with the dichroic 2 A ce n = 437 + 860 setting. This 
mode offers a larger, redder overall wavelength coverage at 
the expense of slightly larger inter-arm gaps (between 5000 
- 6750 A and 8500 - 8750 A). In addition to the UVES data, 
we have obtained ~ 5 hours of Keck HIRES data, which are 
described in more detail by the comprehensive abundance 
analysis of Prochaska et al. (2001). In brief, the C5 decker 



Q2206-199 z-1.9207 N(HI) = 20.6B 



I I 




\ 


J 






iff 


! 





1200 1210 1220 1230 



Q1223+178 z = 2.4661 N(HI)=21.48 




1180 1200 1220 1210 



Wavelength (rest) 

Figure 1. Fits to the damped Lyman a lines for the two QSOs 
studied here. The data are shown in the rest wavelength for the 
absorption rcdshift stated in each panel. 



was implemented providing FWHM w 8 km s _1 resolution 
with the setup covering A « 4750 — 7200 A. The data were 
reduced with the MAKEE software package developed by T. 
Barlow for the extraction of HIRES spectra. The coadded 
spectrum has S/N « 30 per 2 km s _1 pixel. 

To obtain the strictest upper limit on the Co II column 
density for the damped Lyman a system towards Q1223+17, 
we coadded the HIRES and UVES spectra in the region 
spanning the Co II A2012 transition. Each spectrum was 
rebinned to a common wavelength scale and then optimally 
coadded resulting in a final spectrum with S/N » 36 per 2 
km s _1 pixel. 
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4 NEUTRAL HYDROGEN AND COBALT 
COLUMN DENSITIES 

The H I column density of the two QSOs was determined 
by fitting fully damped profiles to the normalised spectra 
using the Starlink package Dipso. The UVES spectrum of 
Q2206-199 was found to have a best fit N(R I) = 4.8±0.3 x 
10 20 cm" 2 and a redshift of z abs = 1.9207, Figure 1. Note 
that this fit is very well constrained by the shape of the 
central trough and damping wings and the iV(H I) is in 
good agreement with previously determined estimates (e.g. 
Pettini et al. 1994) for this DLA. The HIRES spectrum of 
Q1223+17 exhibits a relatively large column density of H I 
which is found to have a best fit of A(H I) = 3±0.3 x 10 21 
cm" 2 and a redshift of z a bs = 2.4661, also in good agreement 
with the Pettini et al. (1994) measurement. 

A detailed analysis of other metal abundances in the 
DLAs towards Q2206-199 and Q1223+17 can be found in 
Prochaska and Wolfe (1997) and Prochaska et al. (2001); 
here we simply concentrate on the detection of Co II. How- 
ever, the analysis of other metal lines plays an important 
role in determining whether lines are blended and distin- 
guishing weak features from the noise. One of the standard 
procedures is to fit Voigt profiles to the absorption system 
which is described as a complex of 'clouds', each defined by 
a Doppler parameter (6- value) , column density (N) and red- 
shift. It has been unanimously found in other high resolution 
DLA studies that singly ionized species such as Fe II, Si II 
and Ni II can be fitted using identical cloud models (e.g. 
Pettini et al 1999). Therefore, line-fitting is often performed 
first for relatively strong lines with several observed transi- 
tions, such as Fe II, in order to fix the b- values and redshifts 
of the constituent components. For the other species, excel- 
lent profile fits can be achieved by allowing only the column 
density to vary as a free parameter, indicating that the singly 
ionized metal lines have similar kinematical structure. Fix- 
ing the cloud model therefore offers a significant advantage 
when fitting weak lines, especially if there is a risk of some 
contaminating absorption from other species at a different 
redshift, e.g. contamination in the Lya forest. 

4.1 Detection Of Cobalt in the DLA Towards 
Q2206 199 

In the case of Q2206— 199, the absorption system is quite 
complex in comparison with other DLA metal line systems. 
The fit was initially constrained using Fe II A2249 and Fe II 
A2260 and refined by performing a simultaneous fit with 
Ni II A1750, Ni II A1741 and Ni II A1709. Together, these five 
transitions can provide excellent constraints on the cloud 
model and, in addition, are two elements that when consid- 
ered along with Co may provide interesting clues to nucle- 
osynthetic processes (see section 5). The b- values and red- 
shifts of each cloud were tied within the fitting procedure 
such that these parameters could vary, but remained iden- 
tical for the two elements. The best fit results for this pro- 
cedure yielded a 6 component cloud model which extends 
over almost 150 km s , see Table 2. Having determined the 
cloud model for Fe II and Ni II, we use this template to fit 
Co II (and the other metal lines), again allowing only the 
column densities to vary. 

The revised /-values of Co (Mullman et al. 1998a,b) 
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Figure 2. Fits to Fe II A2250 and the two Co II lines observed 
towards Q2206— 199 and the non-detection towards Q1223+17. 
In all four panels, the data are shown by the solid histogram and 
the continuum level by the dashed line. For Q2206— 199, the fit 
is shown by the solid curve, and the vertical tick marks in the 
bottom panel mark the expected positions of the two strongest 
components in Q1223+17. The Q2206-199 fit was achieved by 
using Fe II and Ni II transitions to constrain the cloud model 
and then allowing 7V(Co) to vary as a free parameter. The Co II 
A 1466 line is blended in the Lyman a forest and is not used to 
improve the final estimate of ./V(Co). 

show that the two strongest ground state resonance lines 
covered by this spectrum are at rest wavelengths of 1466.21 
A (/=0.031) and 2012.16 A (/=0.037). A third line of sim- 
ilar strength is located at A = 1941.29 A (/=0.034) but 
falls in the gap between the two red arm CCDs. The line 
at 1466 A for a system at z a b s — 1.92 lies in the Lyman a 
forest of the z em = 2.559 QSO. The Co II A2012 line lies 
redward of the QSO's Lyman a emission where the risk of 
mis-identification or contamination is much lower. 

The A2012 line was used to determine the iV(Co) of in- 
dividual components given in Table 2 which also lists the fits 
to the other Fe-peak elements covered by our spectrum. Al- 
though the column density error for any individual compo- 
nent can be quite significant, particularly for closely-spaced, 
blended lines, the error on the total N(X) is usually substan- 
tially smaller. Also, the error clearly depends on the wave- 
length of the transition (since the S/N varies along the spec- 
trum) and is much reduced if several transitions are used in 
the fit. All of these factors are taken into account by VPFIT 
when determining the total column density errors from the 
chi-squared reduction. In addition, we consider the effect of 
the continuum fit on our results, which may be particularly 
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Table 2. Voigt Profile Fit Parameters for Fe-Peak Elements in the DLA Towards Q2206-19. 



5 



Cloud 


Redshift 


b 


Co II 


Ni II 


Logio N(X) 
Fe II 


Zn II 


Cr II 


1 


1.91995 


14.8 


11.54 


13.09 


14.42 


11.88 


12.56 


2 


1.92000 


5.9 


12.51 


13.53 


14.58 


12.10 


12.88 


3 


1.92023 


11.0 


11.81 


13.17 


14.22 


11.90 


12.42 


4 


1.92040 


5.1 


12.05 


12.95 


14.15 


11.87 


12.38 


5 


1.92062 


14.9 


12.45 


13.70 


14.87 


12.45 


13.17 


6 


1.92094 


31.9 


12.61 


13.69 


14.79 


12.41 


13.06 


Total 


system N(X) 


13.09+0.05 


14.23±0.01 


15.36±0.01 


12.95±0.01 


13.63±0.01 



important for Co since a large fraction of the absorption is 
associated with relatively broad, shallow features. We firstly 
note that the continuum fit was done with a low order poly- 
nomial across the flat continuum around the Coll A2012 
line and that any significant departure from this would re- 
quire quite a contrived placement. Secondly, VPFIT may ad- 
just the continuum in order to improve the fit, although no 
change was made in this case. As can be seen in Table 3, the 
relative amounts of Co in the red and blue components are in 
good agreement with the other elements, a sign that at least 
the shape of the continuum is correct. In terms of a system- 
atic offset, in Figure 2 it can be seen that had the continuum 
been placed too low, there would be clear evidence of flux 
above the unity level. If the continuum had been placed too 
high, this would only increase N(Co), thus reinforcing the 
observed overabundance of [Co/Fe]. Nonetheless, we have 
investigated the possible systematic uncertainty induced by 
continuum placement around Co II and other elements and 
determine the effect on derived column densities to be not 
more than 0.10 dex, and typically significantly less than this. 

Over-plotting the fit from Co II A2012 onto the Co II 
A1466 line (Figure 2) shows that although there is evidence 
for weak Co II absorption there is considerable blending, 
probably due to Lyman a forest lines. There is also some 
possibility of blending in the A2012 transition, since it can 
be seen in Figure 2 that an additional component is present 
at the blue edge of the absorption complex that is not fitted 
by the cloud model we adopt (see the top panel of Figure 2 
for a comparison with the Fe II A2250 line) . In order to inves- 
tigate the possible impact of this component on our column 
density determination, we consider the relative contributions 
of the individual clouds to the total absorption seen in the 
other metal lines. Table 3 shows the column density of metal 
enriched gas in the blue (N(X)t) and red (N(X) r ) compo- 
nents, where the blue component is defined to be made up 
of clouds 1 and 2 from Table 2, whereas the red component 
is taken to consist of clouds 3-6. Also listed in Table 3 is 
the fraction of gas in the blue component (/;,) for all of the 
metals measured in this DLA. From these statistics it can 
be seen that, excluding Co, on average 27% of the enriched 
gas resides in the blue component, a value which remains re- 
markably constant for all the ions considered here. For Co, 
however, we see that a slightly higher fraction of gas, 29%, 
resides in the blue component. This small increase may be 
caused by the wing of the contaminating line. We attempted 
to identify the possible cause of contamination, but were un- 
able to associate the extra component with any metal line 
from either of the other 2 DLAs. To show that the effect of 



Table 3. Column Densities in the Blue and Red Components of 
the DLA Towards Q2206-199 



Element, X 


Log N(X) b 


Log N(X) r 


fb 


Co 


12.55 


12.94 


0.29 


Fe 


14.81 


15.22 


0.28 


Ni 


13.67 


14.09 


0.28 


Zn 


12.31 


12.84 


0.23 


Cr 


13.05 


13.50 


0.26 


Si 


15.23 


15.63 


0.28 



the blending is probably negligible, we assume that the frac- 
tional distribution of gas is the same for Co as for the other 
metals and that therefore the AV(Co) = 12.94 represents 73% 
of the total, implying a total N(Co) — 13.08, almost exactly 
what is determined from the original fit. 



4.2 The Upper Limit of Cobalt in the DLA 
Towards Q1223+17 

There is no Co II A2012 absorption visible in the combined 
UVES+HIRES spectrum of Q1223+17, see Figure 2. To 
place an upper limit on the N(Co II) value and in turn the 
Co/Fe ratio, we integrated the apparent column density^ 
over a small region covering the Co II A2012 transition. For 
high S/N data (i.e. in the limit that ln(l — x) « x), this is 
equivalent to deriving the equivalent width over the same 
velocity range and then calculating the column density as- 
suming the linear curve of growth. Although the stronger 
metal line transitions span over 200 km s _1 (e.g. Fe II 
1608), to place the tightest constraints on the Co/Fe ra- 
tio we integrated from —20 km/s to +20 km/s which covers 
the strongest Fe II feature in the velocity profile. Restricting 
our analysis to this velocity range, we find N(Co II) < 12.63, 
N(Fe II)= 15.05 and [Co/Fe] < +0.18 where the limits are 
at the 3cr level. 



* N a (v) = m e CT a {v)/-ne 2 fX where r a (v) = ]n[Ii(v)/I a (v)], It 
and I a are the incident and observed intensity, / is the oscillator 
strength, and A is the rest wavelength of the transition (Savage 
& Sembach 1991). 
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4.3 Corrections Due to Ionization and Dust? 

It is usually assumed in DLA abundance studies that the col- 
umn densities of the singly ionized metal species are a good 
approximation to the total column density of this element, 
based on observations and theoretical modeling (e.g. Viegas 
1995). However, there has been some suggestion (Howk and 
Sembach 1999; Izotov et al. 2000) that ionization corrections 
may in fact be important when deriving DLA abundances. 
Although there is some temptation to invoke ionization cor- 
rections to explain a few unusual abundance trends in DLAs, 
the predictions of these models are often not borne out by 
observations (Pettini et al. 2000; Lattanzio et al 2001). In 
addition, Levshakov et al. (2000) have recently shown that 
the physical conditions adopted by these models are unre- 
alistic when compared with those observed in local DLA 
analogs. Here, we are specifically concerned with the rela- 
tive abundances of Co and Fe. Even in the unlikely event 
that some correction is relevant to the absolute abundances 
as suggested by Izotov et al. (2000), the first and second 
ionization stages of cobalt and iron have very similar po- 
tentials (e.g. 7.88 eV and 7.90 eV for Co I and Fe I) so 
that it is highly unlikely that an ionization correction is re- 
quired when converting N(Co II) and N(Fe II) to [Co/Fe]. 
The same is true for nickel, since the ionization potential for 
Ni I is 7.64 eV. 

The presence of dust in DLAs and its tendency to re- 
move a fraction of the total metals from the gas phase has 
now been well documented (e.g. Pettini et al. 1994; Pettini 
et al. 1997). This renders the task of determining the true 
abundances in these systems somewhat uncertain, although 
several attempts have been made to recover intrinsic abun- 
dances by correcting for Galactic depletion patterns (Vladilo 
1998; Savaglio 2000). In the case of Co, some of the more 
interesting Galactic (stellar) abundance patterns are seen 
not in absolute (e.g. [Co/H]) measurements, but in the val- 
ues of relative abundances like [Co/Fe] and [Co/Cr]. The 
important consideration is therefore the relative depletion 
fractions of these elements, which can be determined from 
local interstellar observations. From Savage and Sembach 
(1996) and Mullman et al. (1998b) we can see that Co, Cr, 
Ni and Fe are all depleted to a very similar level. In addi- 
tion, at relatively low levels of depletion, as is the case for 
Q2206— 199, the small differences between the dust correc- 
tion required for these elements becomes negligible (Pettini 
et al. 2000). 

Therefore, we assume that we can safely investigate rel- 
ative abundance ratios such as [Co/Fe] without the need for 
either corrections for dust or ionization. 



5 DISCUSSION OF COBALT ABUNDANCES 

In Table 4, we list the abundances of the iron peak elements 
(Fe, Ni, Cr, Zn, Co) measured for the two DLAs analysed 
herein and two other systems found in the literature for 
which useful [Co/H] upper limits can be determined. 

5.1 Homogeneity of Stellar Cobalt Abundances 

The cobalt abundances measured in the four DLAs will be 
compared in Section 5.2 with those for stars spanning a 



range of Galactic populations and metallicities (Sneden & 
Gratton 1991; McWilliam & Rich 1994; McWilliam et al. 
1995; Ryan, Norris, & Beers 1996; Prochaska et al. 2000). 
The stellar measurements use solely Co I lines. 

All stellar measurements utilised laboratory gf values 
from Cardon et al. (1982), with the exception of Sneden & 
Gratton who derived solar gf values. These two scales ap- 
pear to be very similar; five lines in common give a mean 
difference log <?/gs90 — log (;/csstw82 = +0.06 dex with a 
standard error in the mean of 0.04 dex. As this is significant 
only at the 1.5<r level, it would be premature to adjust for 
it. Furthermore, by using solar gf values, Sneden and Grat- 
ton have effectively analysed their stars differentially with 
respect to the Sun, so changing the gf values would alter the 
inferred solar abundance as well. We note, nevertheless, that 
even in their most metal-rich stars, [Co/Fe] tends towards 
~ —0.05 at [Fe/H] = 0, which supports the view that their 
Co abundances are marginally low. The problem (assuming 
there is one) may perhaps be traced to differences in the 
assumptions in the solar model used for their solar analy- 
sis and those employed in the other stellar calculations. We 
did, however, correct the three stars of Sneden & Gratton 
affected by the probable 0.47 dex error in the solar log gf 
value of the 4118 A line; see Norris, Ryan & Beers (2001) for 
details. A comparison with the improved gf values of Nitz et 
al. (1999) confirms the accuracy of the Cardon et al. scale: 
log 3/csstw82— log 3/nkwl99 = —0.02 dex and acssTW82 
~ 0.08 dex. We are thus confident that all stars are on a uni- 
form, accurate gf scale. This is quite an achievement since 
the stars span four orders of magnitude in abundance, and 
have been analysed by several independent groups. 

As the DLA Co abundances are derived from ionised 
lines whereas the stellar transitions are neutral, there is (in 
principle) a possibility that the stellar and DLA abundance 
scales differ. However, as the gf scales for both ionisation 
states are tied to modern lifetime measurements, we doubt 
that differences exceed 5-10%, or 0.02 - 0.04 dex. 

Co I lines have extensive hyperfine structure which can 
affect abundance measurements by up to several xO.l dex 
(e.g. Ryan et al. 1996, Fig. 1). The stellar studies cited have 
all accounted for hyperfine structure, which was not always 
the case in earlier analyses. Moreover, work by Pickering 
and colleagues (1998 and priv.comm.) have confirmed that 
the Co II lines we have observed in the DLAs will not be 
broadened sufficiently to affect our results. 

Both Co and Fe are primarily singly ionised in the tem- 
perature range of the stars investigated, so the Co I mea- 
surements are of a minority state. However, because these 
elements have almost identical first ionisation potentials, un- 
certainties in their temperatures or surface gravities will af- 
fect the ionisation degree of Co I just as much as for Fe I. 
Representative errors in [Co/Fe] — see Ryan et al. 1996, 
Table 3 — appear to be ~0.05-0.10 dex (la) for good S/N 
(~100) data, but can be ~ 0.2 dex for poorer data with S/N 
~ 30. It is believed that the lower S/N levels achieved in the 
first studies of the most metal-poor stars may explain most 
of the spread in [Co/Fe] seen in those stars. (See Norris et 
al. 2001 for much improved [Co/Fe] measurements in such 
objects.) 
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Table 4. Iron-peak abundances. Co limits are 3o\ Note that the Ni abundance for HE1 104— 1805 has been increased by a factor of 1.9 
from the published value, based on the revised /-values of Fedchak and Lawler (1999) 



Quasar 


[Fe/H] 


[Ni/H] 


[Co/H] 


[Cr/H] 


[Zn/H] 


Ref 


Q2206-199 


-0.81 


-0.70 


-0.50 


-0.73 


-0.42 


This work 


Q1223+17 


-1.94 


-1.80 


< -1.76 


-1.66 


-1.94 


Prochaska et al. (2001), this work 


Q0302-223 


-1.20 


-1.04 


< -0.87 


-0.96 


-0.56 


Pettini et al. (2000) 


HE1104-1805 


-1.59 


-1.40 


< -1.56 


-1.46 


-1.02 


Lopez et al. (1999) 



5.2 Comparison with Relative Stellar Abundances 

The stellar data discussed in the previous section are plot- 
ted in Figure 3, as a function of [Fe/H]. As discussed in the 
Introduction, Galactic stellar [Co/Fe] ratios show distinct 
— but non-monotonic — trends as a function of metallicity, 
which indicate key differences between the various stellar 
populations. The halo and bulge data, for example, at high 
and low metallicity, share great age and the highest [Co/Fe] 
ratios, whereas later-forming stars in the disk have progres- 
sively lower [Co/Fe] ratios. In addition, the Co detection 
in Q2206— 199 and three upper limits are also plotted, al- 
though for the DLAs [Zn/H] is used as a metallicity indica- 
tor rather than [Fe/H]. The first determination of [Co/Fe] = 
+0.31 ±0.05 in a DLA is consistent with the overabundance 
seen in Galactic halo stars and the transition regime between 
the thick disk and bulge at a metallicity consistent with the 
latter populations. The point plotted for Q2206— 199 is the 
value determined from our Voigt profile fit although, as dis- 
cussed in section 4, this may be a small overestimate if there 
is some blending of the blue components. In addition, Co is 
also overabundant relative to Cr in Q2206— 199, [Co/Cr] = 
+0.23, again consistent with thick disk, moderately metal 
poor halo and metal-rich bulge stars. 

The trend of high [Co/Fe] in metal-poor stars was first 
pointed out by McWilliam et al. (1995) and confirmed by 
Ryan et al. (1996). The observation that stellar [Co/Fe] 
abundances are highest in the oldest populations might sug- 
gest that elevated abundances arise from the fastest evolv- 
ing stars, i.e. the most massive supernova (SN) progenitors. 
However, halo star [Co/Fe] values are not uniformly high, 
and clearly vary as a function of metallicity. McWilliam et 
al. (1995) suspected that the trends were due to a metallicity 
dependence of the yields of elements produced in supernova 
(SN) nucleosynthesis, while Ryan et al. (1996) noted that 
differences in the explosion energy could affect the yields in 
this way. It is unclear whether the metallicity-dependence 
is secondary, being driven by lower Co/Fe yields from later- 
evolving, lower mass SN, or whether it reflects an explicit 
dependence on progenitor metallicity. More than one factor 
may be (and probably is) at work. Umeda et al. (2000) have 
shown that the abundances of these elements do not de- 
pend strongly on the initial stellar metallicity. Nakamura et 
al. (1999) have shown that the abundance trends of Cr, Co 
and Mn observed in the most metal-poor stars can be repro- 
duced when changing the position of the mass cut, but at the 
expense of overproducing Ni, which is only rarely observed 
in the stars (Ryan et al. 1996). More recently, Nakamura 
et al. (2000) have also been able to reproduce the trends, 
though not the absolute values, in Co and Cr by using 
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Figure 3. Relative cobalt abundances as a function of metallicity 
for different Galactic populations and DLAs: open stars — metal- 
poor halo stars from McWilliam et al (1995) and Ryan, Norris, & 
Beers (1996); open circles — moderately metal-poor thin disk and 
halo stars from Sneden & Gratton (1991), with gf changes as de- 
scribed in the text; open triangles — bulge stars from McWilliam 
& Rich (1994); open squares — thick-disk stars from Prochaska 
et al. (2000). The solid circle represents Q2206-199 (with associ- 
ated fitting error) and the down-turned filled triangles are upper 
limits for three other DLAs. The cross in the lower left corner of 
the figure shows the typical error bar for stellar points, whereas 
the error bar in the top right corner shows the maximum effect 
of systematics (continuum placement and N(H I) determination) 
discussed in the text. 



nucleosynthesis models of hypernovae, i.e. supernovae with 
very high explosion energies. These authors remark that as 
well as re-producing Galactic abundances, they can also ex- 
plain the unusual patterns seen in some starburst galaxies, 
such as M82. Indeed, several other observational clues indi- 
cate that star formation in DLAs proceeds in short bursts 
separated by extended quiescent periods, such as low N/O 
and [S,Si/Zn] values (Pettini et al. 2001; Pettini et al 2000) 
and the lack of detected star formation through Ha imaging 
down to levels of ~ 1 Mq yr _1 (Bouche et al. 2000; Kulkarni 
et al. 2001). However, both the models of Nakamura et al. 
(1999; 2000) also predict an overabundance of Ni relative to 
Fe, something that is rarely seen in stars, but quite com- 
monly seen in DLAs (Figure 4). In fact all 4 of the DLAs 
plotted in Figure 3 exhibit [Ni/Fe] ~ 0.2 [j]. However, since at 



t We note that Prochaska & Wolfe (1997) use the old Nill f- 
values, requiring a correction by a factor of 1.9, which brings 
their N(Ni) into very good agreement with ours. 
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least HE1104— 1805 shows no strong overabundance of Co, 
it appears that neither stars nor DLAs exhibit consistent or 
comprehended relative abundances of Co and Ni. Finally, we 
note that since the revision of Ni II /-values, a systematic 
overabundance is now seen for [Ni/Fe] in DLAs, perhaps in- 
dicating that further refinement of the atomic data may still 
be required. 

As detailed supernova nucleosynthesis models have yet 
to provide a convincing explanation of the [Co/Fe] abun- 
dances of stars, we would not presume to be able to solve 
that mystery with the first ever DLA observations. Never- 
theless, by comparing the DLA data with the stellar mea- 
surements, useful insights can be obtained. As noted, the 
[Co/Fe] ratio in the 2 abs = 1-92 DLA of Q2202-199 falls 
near the boundary between bulge and thick disk stars. It 
is seen at a lookback time of ~8.5 Gyr, around the time 
of (or a few Gyr after) the formation of the oldest thick 
disk and bulge objects. However, the possible absence of an 
elevated [Co/Fe] ratio in Q 1223+ 17, with a lookback time 
slightly (~ 0.5 Gyr) longer, emphasises that age alone does 
not explain the elevated [Co/Fe] values. If massive stars do 
produce the elevated [Co/Fe] ratios, then possibly star for- 
mation started later in the z — 1.92 DLA than in the others. 
This postulate could be further explored by considering the 
abundance of Si, an a element. We find [Si/Zn] =—0.07, al- 
though this value is likely be be an underestimate of the 
true Si abundance due to dust depletion. It is possible to 
make a crude adjustment for the depletion fraction by con- 
sidering Zn and Cr ratios, which are found to be solar in 
stars. Therefore, any significant departure from [Zn/Cr] = 
can be explained by the incorporation of Cr into grains. 
For Q2206-199, we determine [Zn/Cr] = 0.35 dex, apply- 
ing this correction to Si gives [Si/Zn] ~ 0.3, although this is 
likely to be an upper limit since Si is less depleted than Cr in 
the local ISM (Savage and Sembach 1996). In any case, it is 
likely that Si exhibits at least a mild overabundance relative 
to Zn, consistent with the expectation if the enrichment of 
this galaxy's ISM has been dominated by the products of 
Type II SN. 

Nakamura et al. (1999, 2000) have achieved partial suc- 
cess in attempting to explain the iron-group element ra- 
tios (Ni excepted) in their nucleosynthesis models, firstly 
by moving the mass-cut and secondly by changing the ex- 
plosion energy. The latter, as they note, partially simu- 
lates the former effect because the zone of complete Si- 
burning is shifted outward in high-energy models compared 
to classical 10 erg models. However, increasing the en- 
ergy has the additional beneficial outcome that it provides 
a solution, for the first time, to the long-standing problem 
of explaining why Ti is overabundant in halo stars. It is 
found that the greater radial extent over which Si-burning 
occurs in the high-energy models necessarily incorporates 
lower density zones, which allows for increased a-rich freeze- 
out-production of 48 Ti. Previous classical SN nucleosynthe- 
sis models always yielded [Ti/Fe] ~ 0. If this finding may 
be considered independent evidence that higher-energy ex- 
plosions are the norm in the more massive SN, i.e. those 
that enriched the extremely metal-poor halo and the 'get- 
rich-quick' bulge, then we should understand the higher 
[Co/Fe] and simultaneously the higher [Ti/Fe] ratios in these 
populations that were enriched by massive-stars. However, 
Prochaska & Wolfe (1997) find that [Ti/Fe] = -0.07 in 
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Figure 4. Ni and Fc abundances in DLAs. Open circles are values 
taken from the literature (Prochaska et al 2001; Pettini et al. 
2000), filled triangles are the DLAs for which a Co upper limit has 
been obtained herein and the filled circle is for the DLA towards 
Q2206— 19 for which Co has been detected. 



Q2206— 199, although considering that Ti is slightly more 
depleted than Fe (Savage and Sembach 1996), this DLA is 
likely to have a roughly solar Ti/Fe ratio. 

Once high-energy SN models are explored more widely 
and included into galactic chemical evolution calculations, 
and a larger range of elements including Co and possibly 
Ti have been measured in more DLAs, the complementary 
constraints presented by DLA and stellar systems can be 
considered jointly. Nevertheless, it is important to appreci- 
ate through all of this that the similarities of DLA [Co/Fe] 
and [Co/Cr] ratios to those of a given Galactic population 
do not necessarily imply that this DLA, or the DLA popula- 
tion in general, represents the precursor to modern day ha- 
los/bulges. Rather, it is indicative of similar physical proper- 
ties such as star formation rate and mix of SN types between 
the populations. 



6 SUMMARY AND CONCLUSIONS 

The distinct trends exhibited by Co in the various Galactic 
stellar populations makes this an interesting element to tar- 
get in DLA abundance studies. However, the weakness of the 
Co II lines renders this a challenging prospect, even for an 
8-m class telescope such as the VLT. Nonetheless, we have 
shown that by selecting DLAs with high column densities of 
Fe towards bright QSOs, the detection of Co is feasible with 
instruments such as UVES and HIRES down to [Co/Fe] ~ 
0. We have achieved the first detection of [Co/H] = —0.50 in 
one of the DLAs studied herein and a sensitive upper limit 
for a second system. In addition, we have determined useful 
upper limits for [Co/H] for two other DLAs from the lit- 
erature. The detection of Co towards Q2206— 199 shows an 
overabundance relative to Fe, [Co/Fe] = +0.31±0.05, at a 
metallicity of approximately \Zq, consistent with both the 
ratios observed in Galactic bulge and thick disk stars. In 
fact, Prochaska et al (2000) noted that there is evidence from 
Galactic bulge/thick disk stellar abundances that these two 
populations may have shared a similar formation history. 
One current limitation of the comparison with bulge stars, 
however, is that the only measurements presently available 
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in the literature were obtained with a 4-m telescope with rel- 
atively low (R ~ 17, 000) resolution (McWilliam and Rich 
1994). Repeating these observations with Keck/HIRES has 
shown that these abundances may undergo significant revi- 
sion, for example [Fe/H] has already been found to be higher 
by ~ 0.2 compared with the 4-m data of lower resolution 
and S/N (Rich and McWilliam 2000). The revised Co abun- 
dances from the HIRES data are not yet available, but will 
eventually provide an interesting comparison for our DLA 
detection. Conversely, the upper limit of [Co/Fe] < +0.18 
that we determine for Q 1223+ 17 is not consistent with the 
bulge star pattern and is more in line with Galactic disk 
abundances, as is the result for HE1104— 1805. 

However, the interpretation of the data presented here 
are limited by statistics and we clearly can not presume 
to understand the star formation history of DLAs based 
on a single detection. In addition, Q2206— 199 is already 
known to exhibit large abundances of other elements which 
are not normally seen in DLAs, such as Ti (Prochaska and 
Wolfe 1997). We can therefore not exclude the possibility 
that Q2206— 199 exhibits peculiar abundances that may not 
be seen in other such systems. Finally, we have deliberately 
chosen DLAs with high N(Fe) for this initial study in order 
to maximise the probability of a Co detection. In order to ex- 
tend this analysis, such a selection bias will eventually have 
to be overcome if we are to truly understand the trends of 
relative abundances. Nonetheless, this is an interesting new 
result and the study of Co in DLAs clearly invites further 
study. 

Our understanding of the processes that drive abun- 
dance trends in stars has only truly become possible thanks 
to large datasets, since intrinsic scatter and errors will dom- 
inate small samples. Therefore, more DLA data points are 
clearly required in order to draw firm conclusions about the 
star for mation histories implied by the observed abundance 



trends. fPerhaps most importantly wp have shown that, with 

current 8-m class telescopes and high resolution echelle spec- 
trographs that such measurements are not only possible, but 
also that Co may be a rewarding element to study in the 
context of testing nuclesynthetic models and understanding 
galactic star formation histories. 
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